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Plane inertial waves are generated using a wavemaker, made of oscillating stacked plates, in 
a rotating water tank. Using particle image velocimetry, we observe that, after a transient, the 
primary plane wave is subject to a subharmonic instability and excites two secondary plane waves. 
The measured frequencies and wavevectors of these secondary waves are in quantitative agreement 
with the predictions of the triadic resonance mechanism. The secondary wavevectors are found 
systematically more normal to the rotation axis than the primary wavevector: this feature illustrates 
the basic mechanism at the origin of the energy transfers towards slow, quasi two-dimensional, 
O ' motions in rotating turbulence. 



o 



>> 
i 

q 

o 

c/a 
>> 

Oh. 




2 



I. INTRODUCTION 

Rotating and stratified fluids support the existence of two classes of anisotropic dispersive waves, called respectively 
inertial and internal waves, which play a major role in the dynamics of astrophysical and geophysical flows.— ~— These 
waves share a number of similar properties, such as a group velocity normal to the phase velocity. Remarkably, in 
both cases, the frequency of the wave selects only its direction of propagation, whereas the wavelength is selected by 
other physical properties of the system, such as the boundary conditions or the viscosity.— »^ 

Most of the previous laboratory experiments on inertial waves in rotating fluids have focused on inertial modes or 
wave attractors in closed containers^— whereas less attention has been paid to propagative inertial wave beams. 
Inertial modes and attractors are generated either from a disturbance of significant size compared to the container^ 
or more classically from global forcing.— ~— Inertial modes are also detected in the ensemble average of rotating 
turbulence experiments in closed containers ! 13 ! 14 On the other hand, localized propagative inertial wave beams have 
been investigated recently in experiments using particle image velocimetry (PIV) i 15 ' 16 

A monochromatic internal or inertial wave of finite amplitude may become unstable with respect to a parametric 
subharmonic instability— ~— This instability originates from a nonlinear resonant interaction of three waves, and 
induces an energy transfer from the primary wave towards two secondary waves of lower frequencies. This instability 
has received considerable interest in the case of internal gravity waves^ because it is believed to provide an efficient 
mechanism of dissipation in the oceans, by allowing a transfer of energy from the large to the small scalesi^— 

Parametric instability is a generic mechanism expected for any forced oscillator. A pendulum forced at twice its 
natural frequency provides a classical illustration of this mechanism. Here, the "parameter" is the natural frequency 
of the pendulum, which is modulated in time through variations of the gravity or pendulum length. Weakly nonlinear 
theory shows that the energy of the excitation, at frequency <7q, is transferred to the pendulum at its natural frequency 
cto/2, resulting in an exponential growth of the oscillation amplitude. 

In the case of inertial (resp. internal) waves, the "parameter" is now the so-called Coriolis frequency / = 251, 
with f2 the rotation rate (resp. the Brunt- VaisaTa frequency ./V). In the presence of a primary wave of frequency <7o, 
this "parameter" becomes locally modulated in time at frequency erg, and is hence able to excite secondary waves 
of lower natural frequency. However, here a continuum of frequencies can be excited, so that the frequencies o\ and 
01 of the secondary waves are not necessarily half the excitation frequency, but they nevertheless have to satisfy the 
resonant condition g\ + a-2 = ctq. Interestingly, in the absence of dissipation, the standard pendulum- like resonance 
o"i = (72 = co/2 is recovered both for inertial and internal waves, and the corresponding secondary waves have 
vanishing wavelengths. 20 Viscosity is responsible here for the lift of degeneracy, by selecting a maximum growth rate 
corresponding to finite wavelengths, with frequencies a\ and a 2 splitted on both sides of o~o/2£^ 

The parametric subharmonic instability has been investigated in detail for internal gravity wavesj 20 ' 24 On the other 
hand, this instability mechanism has received less attention in the case of pure inertial waves (i.e., in absence of 
stratification) , probably because of the lower importance of rotation effects compared to stratification effects in most 
geophysical flows. It has been observed in numerical simulations of inertial modes in a periodically compressed rotating 
cylinder! 10 ' 11 To our knowledge, parametric instability in the simpler geometry of plane inertial waves has not been 
investigated so far, and is the subject of this paper. A fundamental motivation for this work is the key role played by 
triadic interactions of inertial waves in the problem of the generation of slow quasi-2D flows in rotating turbulence 
The parametric subharmonic instability indeed provides a simple but nontrivial mechanism for anisotropic energy 
transfers from modes of arbitrary wavevectors towards lower frequency modes, of wavevector closer to the plane normal 
to the rotation axis (i.e., more "horizontal" by convention). Note that this nonlinear mechanism may however be in 
competition with a linear mechanism — the radiation of inertial waves along the rotation axis — which has also been 
shown to support the formation of vertical columnar structures^ The relative importance of these two mechanisms 
is governed by the Rossby number, defined as Ro = (T n if7) _1 , with fi" 1 the linear timescale and r n i = L/U the 
nonlinear timescale based on the characteristic velocity U and length scale L. In rotating turbulence with Ro <C 1, 
the anisotropy growth should hence be dominated by the nonlinear triadic interactions, whereas for Ro = 0(1) both 
mechanisms should be at play. 

In this paper, we report the first experimental observation of the dcstabilization of a primary plane inertial wave 
and the subsequent excitation of subharmonic secondary waves. To produce a plane inertial wave of sufficient spatial 
extent, and hence of well-defined wavevector ko, we have made use of a wave generator already developed for internal 
waves in stratified fluidsi22r— Wave beams of tunable shape and orientation can be generated with this wavemaker. We 
show that, after a transient, the excited plane wave undergoes a parametric subharmonic instability. This instability 
leads to the excitation of two secondary plane waves, with wavevectors which are systematically more "horizontal" 
than the primary wavevector. We show that the predictions from the resonant triadic interaction theory for inertial 
waves, as described by Smith and Walcffe^ are in excellent agreement with our experimental results. In particular, 
the frequencies and wavenumbers of the secondary waves accurately match the expected theoretical values. 
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FIG. 1. (Color online) Schematic representation of the wave generator. The excited plane inertial wave has a frequency <ro, a 
downward phase velocity, a negative helicity (so = —1), and propagates at an angle 9 = cos - (cro//), with / = 2f2 the Coriolis 
parameter. 



II. INERTIAL PLANE WAVE GENERATION 



A. Structure of a plane inertial wave 



We first briefly recall the main properties of inertial waves in a homogeneous fluid rotating at a constant rate Q. 
In the rotating frame, the restoring nature of the Coriolis force is responsible for the propagation of the inertial 
waves, for frequencies a < /, where / = 2£l is the Coriolis parameter. Fluid particles excited at frequency a describe 
anticyclonic circles in a plane tilted at an angle 9 = cos" 1 (u//) with respect to the horizontal, and the phase of this 
circular motion propagates perpendicularly to this tilted plane. 

The equations of motion for a viscous fluid in a frame rotating at a rate O = //2 around the axis z are 

d t u + (u- V)u = — Vf>- fe z x u + tA7 2 u, (1) 
P 

V-u = 0, (2) 

where u = (u x ,u y ,u z ) is the velocity field in cartesian coordinates x = (x, y, z). In the following, we restrict to the 
case of a flow invariant along the horizontal direction y. The fluid being incompressible, the motion in the vertical 
plane (x, z) may be described by a streamfunction ip{x, z), such that u = (d z ip^u y , —d x ip). Neglecting viscosity, the 
linearized equations for small velocity disturbances are 

d t d z ip = --d x p + fu y , (3) 

d t u y = -fd z tp, (4) 

-d t d x ip = — d z p. (5) 
P 

These equations may be combined to obtain the equation of propagation for inertial waves, 

d tt (d xx + d zz )i> + f 2 d zz i> = 0. (6) 

Considering a plane wave solution of frequency a and wavevector k = (k, 0, m), 

^(x,z,i)=V'oe l(k - x ^ t) +c.c. (7) 

(where c.c. means complex conjugate), we obtain the anisotropic dispersion relation for inertial waves 

771 

a = sf — = s/cos6>, (8) 

with k — (k 2 + to 2 ) 1 / 2 , s = ±1, and 9 the angle between k and the rotation axis (see Fig. [T]). We see from Eq. © 
that a given frequency a lower than / selects a propagation angle ±0, without specifying the norm of the wavevector 
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k. The corresponding velocity field is given by 

im^e' (fa+m2 " fft) +c.c. (9) 

sk%1)q e i{kx+mz - at) +c.c. (10) 

-ikilj e i( - kx+mz - at '> +c.c. (11) 

We recover here that the fluid particles describe anticyclonic circular motions in tilted planes perpendicular to k, as 
sketched in Fig. [T] The wave travels with a phase velocity c v = crk/K 2 and a group velocity c g = Vk<7 normal to c<~. 
The vorticity u) = V x u, given by 

u> = —sku, (12) 

is associated to the shearing motion between planes of constant phase. Because the velocity and vorticity are aligned, 
inertial waves are also called helical waves, and the sign s in Eq. (|5J) identifies to the sign of the wave helicity u • u>, 
with s = +1 for a right-handed wave and s = — 1 for a left-handed wave. For instance, in the classical St. Andrew's 
wave pattern emitted from a point source^ the two upper beams are right-handed and the two lower beams are 
left-handed, although the fluid motion is always anticyclonic. 



u x = 

Uy = 

u z = 



B. Generation of a plane inertial wave 



In order to generate a plane inertial wave, we have made use of a wavemaker, introduced by Gostiaux et al.^ 
which was originally designed to generate internal gravity waves (see Mercier et al£^ for a detailed characterization 
of the wavemaker). This wavemaker consists in a series of oscillating stacked plates, designed to reproduce the fluid 
motion in the bulk of an internal gravity wave invariant along y. The use of this internal wave generator for the 
generation of inertial waves is motivated by the similarity of the spatial structure of the two types of waves in the 
vertical plane (x, z). However, the fluid motion in the internal wave is a simple oscillating translation in the direction 
of the group velocity, whereas fluid particles describe anticyclonic circular translation in the case of inertial waves. As 
a consequence, the oscillating plates of the wavemaker only force the longitudinal component of the circular motion 
of the inertial waves, whereas the lateral component is let to freely adjust according to the spatial structure of the 
wave solution. 

The wavemaker is made of a series of 48 parallelcpipedic plates stacked around a helical camshaft, with the ap- 
propriate shifts between successive cames in order to form a sinusoidal profile at the surface of the generator. We 
introduce the local coordinate system (£, y, rj), tilted at an angle 9 about y, where £ is along the wave propagation and 
X] is parallel to the camshaft axis (see Fig.[T]). The group velocity and the phase velocity of the wave are oriented along 
£ and r\ respectively. As the camshaft rotates at frequency <t , the plates, which are constrained in the y direction, 
oscillate back and forth along £. The sign of the rotation of the helical camshaft selects the helicity of the excited 
wave, and hence an upward or downward phase velocity. In the present experiment, the rotation of the camshaft is 
set to produce a downward phase velocity, resulting in a left-handed inertial wave of negative helicity sq = — 1. 

The cames are 14 cm wide in the y direction, and their eccentricities are chosen to produce a sinusoidal displacement 
profile, £o(?7) = £o sin(Ko^); of wavelength A = 2tt/ko = 7.6 cm and amplitude £ = 0.5 cm at the center of the beam. 
The wave beam has a width 30.5 cm with a smooth decrease to at the borders, and contains approximately 4 
wavelengths. The generator is only forcing the £ component of the inertial wave, and the y component is found to 
adjust according to the inertial wave structure after a distance of order of 2 cm. 

The wavemaker is placed in a tank of 120 cm length, 80 cm width and 70 cm depth which is filled with 58 cm of 
water. The tank is mounted on the precision rotating platform "Gyroflow" of 2 m in diameter. The angular velocity 
f2 of the platform is set in the range 1.05 to 3.15 rad s _1 , with relative fluctuations less than 10~ 3 . A cover is 

placed at the free surface, preventing from disturbances due to residual surface waves. The rotation of the fluid is set 
long before each experiment (at least 1 hour) in order to avoid transient spin-up recirculations and to achieve a clean 
solid body rotation. 

The propagation angle 9 of the inertial wave is varied by changing the rotation rate of the platform, while keeping 
the wavemaker frequency constant , <jq = 1.05 rad s^ 1 . This allows to have a fixed wave amplitude oq £ q = 0.52 cm s _1 
for all angles. The Coriolis parameter has been varied in the range / = 1.004 Co to 3oo, corresponding to angles 9 
from 5° to 70°. For each value of the rotation rate, the axis of the wavemaker camshaft is tilted to the corresponding 
angle 9 = cos _1 (cro//), in order to keep the plate oscillation aligned with the fluid motion in the excited wave. As 
a consequence, the efficiency of the forcing should not depend significantly on the angle 9. For each experiment, the 
fluid is first reset to a solid body rotation before the wavemaker is started. 
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(a) t = 2T (b) t = 7T 
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FIG. 2. (Color) Horizontal velocity field after 2 and 7 periods from the start of the wavemaker for (Jo/ f = 0.84. The wavemaker 
is on the top-right, forcing a wave propagating along c g with a phase propagating along c v . 

C. PIV measurements 

Velocity fields are measured using a 2D particle image velocimetry (PIV) syste m 32 ' 33 mounted on the rotating 
platform. The flow is seeded by 10 /jm tracer particles, and illuminated by a vertical laser sheet, generated by a 
140 mJ Nd:YAG pulsed laser. A vertical 59x59 cm 2 field of view is acquired by a 14 bits 2048x2048 pixels camera 
synchronized with the laser pulses. For each rotation rate, a set of 3200 images is recorded, at a frequency of 4 Hz, 
representing 24 images per wavemaker period. This frame rate is set to achieve a typical particle displacement of 
5 to 10 pixels between each frame, ensuring an optimal signal-to-noise ratio for the velocity measurement. PIV 
computations are performed over successive images, on 32x32 pixels interrogation windows with 50% overlap. The 
spatial resolution is 4.6 mm, which represents 17 points per wavelength of the inertial wave. 

Figure [5] shows typical instantaneous horizontal velocity fields after 2 and 7 periods T = 27t/cto from the start of 
generator, for an experiment performed with <7o// = 0.84. A well defined truncated plane wave propagates downward, 
making an angle 9 = cos~ 1 (cro//) — 34° to the horizontal. The front of the plane wave is propagating at a velocity 
8.3±0.6 mm s , which agrees well with the expected group velocity c g = fshx6/K = 8.5 mm s _1 . The phase velocity 
is downward, normal to the group velocity, and also agrees with the expected value c v = <jq/k = 12.7 mm s" 1 . 

Two sources of noise have been identified, which can be seen in the temporal energy spectrum of the velocity 
fields (Fig. [31 described in the next subsection): an oscillatory motion at frequency a = Q = 0.5/, due to a residual 
modulation of the rotation rate of the platform, and slowly drifting thermal convection structures at frequency 
er — > 0, due to slight temperature inhomogeneities in the tank. Both effects contribute to a velocity noise of order of 
0.2 mm s _1 , i.e. 25 times lower than the wave amplitude close to the wavemaker. This noise could be safely removed 
using a temporal Fourier filtering of the velocity fields at the forcing frequency a®. This filtering however fails in the 
particular case where ctq — Q, for which the mechanical noise of the platform cannot be filtered out of the inertial 
wave signal. 

The wavemaker is found to successfully generate well defined plane waves for frequencies oo > 0.65/. For lower 
frequency, i.e. for steeper angle of propagation [0 = cos -1 (cro//) > 50°], the wave pattern shows significant departure 
from the expected plane wave profile, which may be attributed to the interference of the incident wave with the 
reflected wave on the bottom of the tank. 



III. SUBHARMONIC INSTABILITY 

A. Experimental observations 

After a few excitation periods, the front of the inertial wave has travelled outside the region of interest, and the 
inertial wave can be considered locally in a stationary regime. However, after typically 15 wavemaker periods (the 
exact value depends on the ratio (Jo/ f), the inertial wave becomes unstable and show slow disturbances of scale slightly 
smaller than the excited wavelength. 

We have characterized this instability using Fourier analysis of the PIV time series. We compute, at each location 
(x, z) of the PIV field, the temporal Fourier transform of the two velocity components over a temporal window Ai, 

! rto+At 

u a (x,z) = -= u(x,z,t)e iat dt. (13) 

V 27T J t 
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FIG. 3. (Color online) Temporal energy spectra for two experiments performed at rotation rate f2 = 0.63 rad s _1 with 
(continuous line) and without (dashed line) the wavemaker operating at (To// = 0.84. The spectrum with the generator 
working has been computed on the time interval between 24 and 116 periods after the start of the generator. The peak at 
a/f = 0.5 present in the two spectra is the trace of the mechanical noise of the platform at the rotation frequency a = fi, 
whereas the low frequencies are due to thermal convection effects (see text). 
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TABLE I. Frequencies of the secondary waves ai/f and CT2//, determined from the peaks in the temporal energy spectra, as a 
function of the frequency of the primary wave (Jo/ S- The uncertainty for <J\j j and 02// is ±0.03. 



The temporal energy spectrum is then defined as 

E{a) = (|u.| 2 k z , (14) 

where (-) x ,z is the spatial average over the PIV field. 

If we compute E(a) over a temporal window At spanning a few excitation periods, we observe, as to is increased, 
the emergence of two broad peaks at frequencies smaller than the excitation frequency cto, suggesting the growth of a 
subharmonic instability. These two subharmonic peaks can be seen in Fig. [3l for an experiment performed at rotation 
rate ft = 0.63 rad s _1 with the wavemaker operating at (To/ f = 0.84. Here, the temporal window At is chosen equal 
to 92 wavemaker periods, yielding a spectral resolution of Act = 2w/At ~ 9 x 10~ 3 /. The two secondary peaks 
are centered on cti// = 0.25 ± 0.03 and 02/ / = 0.59 ± 0.03, and their sum matches well with the forcing frequency 
Co// = 0.84, as expected for a subharmonic resonance. The significant width of the secondary peaks, of order 0.07/, 
indicates that this resonance is weakly selective. This broad-band selection will be further discussed in Sec. IIV B[ 

The subharmonic instability of the primary wave is found for all forcing frequencies cto ranging from 0.65/ to /; the 
measured frequencies a\^ are given in Tab. HI The absence of clear subharmonic instability at lower forcing frequency 
may be due to an intrinsic stability of the primary wave for cto < 0.65/, or to the low quality of the plane wave at 
steep angles because of the interference with the reflected wave beam on the bottom of the tank. 

Using temporal Hilbert filtering ; 30 ' 34 the spatial structure of the wave amplitude u D (x) and phase c/?(x, t) = k x at 
can be extracted for each secondary wave. The procedure consists in (i) computing the Fourier transform u CT (x, z) of 
the velocity field according to Eq. (fl"3|) . with a temporal window At of at least 42 excitation periods; (ii) band-pass 
filtering u CT (a;, z) around the frequency of interest <j\ or 02 with a bandwidth of 5a = 2.0 10 -2 /, but without including 
the associated negative frequency; (iii) reconstructing the complex velocity field by computing the inverse Fourier 
transform (including a factor 2, which accounts for the redundant negative frequency, in order to conserve energy), 



u ff (x,t)=u (x) 



(15) 
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FIG. 4. (Color) Hilbert filtered vertical velocity (a) and phase (b) of the primary wave at o"o// = 0.84, and phase of the Hilbert 
filtered first [(c), ai/f — 0.25] and second [(d), 0^2/ f = 0.59] subharmonic waves. The phase is displayed only where the wave 
amplitude is larger than 1.3 10 _1 <ro£o for (b) and 7.7 10 -3 cto£o f° r ( c ) an d (d). In (a), the square in dashed lines indicates the 
region where the primary wave amplitude Aq has been measured. 



The physical velocity field is finally given by Re(uj})- The wave amplitude u D and phase field <p are finally obtained 
from the Hilbcrt-filtercd field uh- 

In Figs. |4|c) and (d), for the experiment at ctq// = 0.84, we show the maps of the phase of the secondary waves, 
extracted from Hilbert filtering at frequencies o\ and oi respectively. It is worth to note, as can be verified from 
Fig. [3l that the corresponding typical velocity amplitude is at least ten times smaller than for the primary wave [see 
Fig. 5[a)]. The spatial structures of the phase of these secondary waves are not as clearly defined as for the primary 
wave [Fig. HJb)]. In particular, dislocations can be distinguished in the phase field. The finite extent of the primary 
wave and its spatial decay due to viscous attenuation are probably responsible for this departure of the secondary 
waves from pure plane waves. It is also important to note that the monochromaticity of the first subharmonic wave 
[Fig. Etc)] is affected by interferences with its reflection on the wavemaker which is due to the fact this secondary 
wave is propagating toward the wavemaker. However, to a reasonable degree of accuracy, the two secondary waves 
can be considered locally as plane waves, characterized by local wavevectors ki and k2. 



B. Helical modes 

The approximate plane wave structure of the two secondary waves suggests to analyze the instability in terms of 
a triadic resonance between the primary wave, of wavevector ko, and the two secondary waves, of wavevectors ki 
and k2. This resonance may be conveniently analyzed in the framework of the helical decomposition, introduced by 
Walcffc , 35 ' 36 which we briefly recall here. 

Helical modes have been introduced as a general spectral decomposition basis, which is useful to analyze the energy 
transfers via triadic interactions. Although this decomposition also applies for non-rotating flows, it is particularly 
relevant for rotating flows, because incrtial plane waves have exactly the structure of helical modest Any velocity 
field can actually be decomposed as a superposition of helical modes of amplitudes A Sk (k, t), 

u( x ,t)=J2 E ^(k^h.JkK^-^, (16) 

k s k =±l 

where <r k k is the frequency associated to a plane wave of wavevector k and helicity sign Sk- The helical mode h Sk (k) 
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is normal to k (by incomprcssibility) , and given by 

. . k k x k x / , 

h - (k) = R x M^ + JSk M^' (17) 

where Sk = ±1 is the sign of the mode helicity^i Injecting the decomposition (fTBf into the Navier-Stokes equation (TTJ) 
yields 

d 



Ak = 1 C^^J^A^e i ^ +a ' +a ^ t , (18) 



with stars denoting complex conjugate, and A^, Uk being short-hands for A Sk (k,t), <r k . In Eq. (p~8|) . the sum is to be 
understood over all wavevectors p and q such that k + p + q = and all corresponding helicity signs s p and s q . In 
the following, the equation k + p + q = will be referred to as the spatial resonance condition for a triad of helical 
modes. The interaction coefficient is given by 



= I \ SaKa - suk p ] (h: p ( P ) x h: q ( q )) • h: k (k). (19) 



kpq ~ 2 1 q q P 



C. Resonant triads 



The helical mode decomposition (|16|) applies for any velocity field, containing an arbitrary spectrum of wavevectors. 
We restrict in the following the analysis to a set of three interacting inertial waves of wavevectors (k, p,q). Equa- 
tion (fT5|) shows that the amplitude of the mode of wavevector k is related to the two other modes p and q according 
to 

+ vhA A k = CkA^A^ e H°*+<r*+"*)* (20) 



where Ck is short-hand for C^ p s * s ' i = C'^^ Sp . Cyclic permutation of k, p and q in Eq. (|2U)) gives the two other 
relevant interaction equations between the three waves. We further restrict the analysis to plane inertial waves 
invariant along y (i.e., k • e y = 0). The three considered helical modes (fT7|) therefore reduce to 

h Sr (r) = is r e y , (21) 



where r stands for k, p or q. From Eq. (|2Tj) . the interaction coefficients (1191) can be explicitly computed, 

C k = — [iTlpkq - m q fc p ] [K q - Kp + SqS k KqK k - S p S k K p K k ], (22) 

and similarly for the two cyclic permutations. 

Since in Eq. (f2"U)l and in its two cyclic permutations the A T (t) coefficients have to be understood as complex velocity 
amplitudes evolving slowly compared to wave periods 2-7r/<7 r , temporal resonance is needed in addition to spatial 
resonance for the left-hand coefficients A T to be nonzero. Using 0, 1, 2 for reindexing the three waves k, p and q, this 
leads to the triadic resonance conditions 

k + k x + k 2 = 0, (23) 

cr + (Jl + cr 2 = 0. (24) 

We consider in the following that only the primary wave Aq, of given frequency do, wavevector ko = (fco,mo) and 
helicity sign sq, is present initially in the system (i.e., ^1.2(0) = 0). The two secondary waves (si, a±, ki) and (S2, 
<7 2 , k 2 ) which could form a resonant triad with the primary wave may be determined using the resonance conditions 
(|2"5|) and (|2"4"|) . From the dispersion relation for inertial waves ((S]), the resonance conditions lead to 

m mi m + mi 

Sp h .Si = — .s 2 — 0. (25) 

Vfco+ m o V k i+ m l V(fc + A;i) 2 + (m + mi) 2 

For a given primary wave (so, ko, mo), the solution of this equation for each sign combination (so, si, s 2 ) is a curve in 
the (fci,mi) plane (see Fig. [5|). Without loss of generality, once we have taken so = —1 (which corresponds to the 
experimental configuration), it is necessary to consider four sign combinations: ( — ,—,—), (—,+,—), (—,—,+) and 
( — , +, +). Notice that the three first combinations always admit solutions, whereas the fourth one, (— , +, +), admits 
a solution only if |m | < Ko/2, i-C. > 60°. The exchange of k x and k 2 keeps the (— , — , — ) and ( — , +, +) resonances 
unchanged, but exchanges the (— , — , +) and ( — , +, — ) resonances. Eventually, three independent sign combinations 
remain: (-,-,-), (— ,=F,±) and (-,+,+). 
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FIG. 5. (Color online) Resonance curves for the primary waves (a) [so = — ljffo = 0.84/, kq = 0.82 rad cm -1 ] and (b) 
[so = —l,cro = 0.99/, fto = 0.82 rad cm - ]. The curves represent the location of ko + ki = (ko + kx,mo + mi) satisfying 
Eq. (|25l) for the 3 possible combinations of signs. The wavevectors measured experimentally are shown using arrows. The circle 
is the theoretical prediction for the location of ko + ki obtained from the maximum growth rate criterion, determined using 
the experimental primary wave amplitude [Ao = 0.29 ±0.07 cm s~ for (a) and Ao = 0.34 ±0.11 cm s" 1 for (b)]. The diameter 
of the circle measures the uncertainty of the prediction due to the uncertainty on the wave amplitude Aq. 



D. Experimental verification of the resonance condition 



The predictions of the triadic resonance theory are compared here with the measured wavevectors of the secondary 
waves. Figure [5] shows the theoretical resonance curves for two forcing frequencies, &o/f = 0.84 and 0.99. For both 
curves, helicity sign and wavenumbcr of the primary wave are chosen according to the experimental values, so = ~~ 1 
and Ko = 0.82 rad cm -1 . 

For both frequencies ao considered here, only the three first sign combinations admit solutions. The (—,—,—) 
combination gives a closed loop, whereas the two others, (— ,±,±), give infinite branches, tending asymptotically 
to constant angles. The limit of large secondary wavevectors is such that |cti| = |cr 2 1 = |<7o|/2: when a wave ko 
excites two waves of wavelength A <C 2tt/kq, both secondary waves have frequency <7o/2, with opposite wavevectors, 
leading to a stationary wave pattern. However, such large wavenumbers are prevented by viscosity, as will be shown 
in Sec. ITVAl 

Figure [5] also shows the measured secondary wavevectors ki and k2. These wavevectors are obtained from the phase 
fields ipi t 2 extracted by Hilbcrt filtering, using 

ki,2 = V^i, 2 . (26) 

These measurements are then averaged over regions of about (130 mm) 2 where the secondary waves can be considered 
as reasonably spatially monochromatic. It must be noted that a same plane wave can be equivalcntly described by (s, 
a > 0, k) and (s, —a < 0, — k). Since we always consider primary waves with positive frequency ao > 0, according to 
Eq. (|24|) . the subharmonic frequencies a\p have to be taken negative. As a consequence, the Hilbert filtering should 
be performed for the negative peaks in the temporal Fourier transform, in order to produce phase fields with the 
appropriate sign. Practically, the Hilbcrt filtering has been performed around the positive peaks — eri.2, and the signs 
of the measured wavevectors have been changed accordingly. 

The secondary wavevectors ki = (fci,mi) and k2 = (£2,7712) measured experimentally, shown in Fig. are in 
good agreement with the triadic condition (f23|) . forming a triangle such that ko ± ki + k2 = 0. Moreover, the apex 
of the triangle, at ko ± ki, falls onto one of the three resonant curves. The selected resonant curve corresponds 
to the sign combination (—,+,—), in agreement with the observed experimental hclicities. We actually verify that 
Si = a\Ki/ /mi is positive (o\ < and mi < 0) and that S2 = a-i^ij jfni is negative (02 < and m-i > 0), confirming 
the (—,+)—) nature of the experimental resonance. 

Interestingly, the shape of the triangle k ± ki ± k 2 = in Fig. [5] indicates that the group velocity of the secondary 
wave ki is oriented towards the wavemaker. Indeed, we recall that, for a given wavevector k, the group velocity c g 
is normal to k, and the vertical projections of c g and k are oriented in the same direction if a > and in opposite 
directions if a < 0. Accordingly, Fig. [5] shows that c g o and c g 2 are oriented downward, pointing from the wavemaker 
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towards the bottom of the tank, whereas c g i is oriented upward, pointing towards the wavemaker. As a consequence, 
the secondary wave k! is fed by the primary wave, but releases its energy back to the wavemaker. 

For all the primary wave angles for which the instability is observed, the secondary waves are systematically such 
that |<7i| and \a^\ are lower than | cr 1 . The dispersion relation hence yields secondary wavevectors ki,2 more horizontal 
than ko , as illustrated in Fig. [5] This property, which actually follows from the conservation of energy and helicity^ 
illustrates the natural tendency of rotating flows to transfer energy towards slow quasi-two-dimensional modes. If 
the process is repeated, as in rotating turbulence, the energy becomes eventually concentrated on nearly horizontal 
wavevectors, corresponding to a quasi-2D flow, with weak dependence along the rotation axisi 26 ' 27 



IV. SELECTION OF THE MOST UNSTABLE RESONANT TRIAD 



A. Maximum growth rate criterion 

In order to univocally predict the resonant secondary waves, a supplementary condition must be added to Eq. (|25|) : 
we assume that the selected resonant triad is the one with the largest growth rate. Going back to the wave interaction 
equations (|20|) associated to the temporal resonance condition (|24)l . the amplitudes of the secondary waves are governed 
by 

- 1 ± = C 1 AIAI-vk\A u (27) 
^=C 2 A* A$-vk 2 2 A 2 , (28) 



with C\ t 2 given by Eq. (|22l) taking k = ki ; 2 (see also Appendix A in Ref. 1251 ). Solving this system with initial 
conditions ^1,2(0) = 0, and assuming that A remains almost constant at short time, lead to the solutions 



A h2 {t)=B l>2 {e^ t -ei- t ), (29) 



where the growth rates j± write 



In the following, we consider the primary wave amplitude as real without loss of generality, so \Aq\ = Aq. 

The coefficient 7_ is always negative, so the stability of the system is governed by the sign of 7+, which we simply 
note 7 in the following. Interestingly, this growth rate 7 depends on the amplitude Aq of the primary wave. As a 
consequence, the primary wave is unstable with respect to a given set of secondary waves, selected by the resonance 
condition and unequivocally denoted by K\, only if A exceeds the threshold A c (ki) = vk 1 k 2 / 'v r C\C 2 in which case 
7(^1) > 0. In other words, for a given couple of secondary waves (denoted by n{) to be possibly growing, the Reynolds 
number based on the primary wave, Reo = Aq/(kqv), must exceed a critical value Re c (ni) = A c {k\) / {kov) for the 
onset of the parametric instability. This critical Reynolds number is actually an increasing function of K\ and tends 
to zero as Ki — > 0, showing that whatever the value of -Reo, there is always a continuum of resonant triads with 
Reo > Re c (ni), i.e. with a positive growth rate. The main consequence is that, whatever the value of Reo-, the most 
unstable triad always has a positive (maximum) growth rate and the parametric instability does not have any Reo 
threshold to proceed. 

If viscosity can be neglected, Eq. (|30|) reduces to 7 = \/ C\C 2 Ao- In the limit of large secondary wavenumbcrs 
Ki 2 3> K 0i one has ki ~ — k2, and the growth rate 7 is found to tend asymptotically toward a maximum valued 
i.e., the selected secondary waves have frequency exactly half the forcing frequency. Taking viscosity into account 
reduces the growth rate of the large wavenumbcrs, and hence selects finite wavenumbcrs. Equation (|30[) indicates 
that larger wavenumbers are selected for larger primary wave amplitudes A Q and/or lower viscosity, i.e. for larger 
Reynolds number Reo. 



B. Selection of the most unstable wavenumbers 



In Fig. |6l the predicted growth rates 7 are plotted for the three possible sign combinations, for the primary wave 
defined by s = — 1, (To = 0.84/, kq = 0.82 rad cm" 1 . These growth rates have been computed using the primary wave 
amplitude averaged over the area where the secondary wavevectors have been measured (see the square in Fig. EJa)), 
Aq = 0.29 cm s _1 . For the 3 types of resonance, the growth rates tend to zero when fci — » — ko/2 and k\ —> 00 
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FIG. 6. (Color online) Growth rates 7 as a function of ki, computed from Eq. (|30[). for the three possible resonances for a 
primary wave (so = — 1,<to = 0.84/, Ko = 0.82 rad cm -1 ). The growth rates have been computed using the average value 
Ao = 0.29 cm s -1 for the primary wave amplitude. For resonance (—,+,—), an additional curve (continuous line) has been 
computed using a wave amplitude 25% larger. 



(because of viscosity) . If the secondary waves ki and 1*2 are exchanged, which amounts to exchange the ( — ,—,+) 
and (— , +, — ) resonances, the same growth rates are obtained: the curves for (— , — , +) and (— , +, — ) arc symmetrical 
with respect to ko/2. 

Interestingly, the growth rate is positive for a broad range of wavenumbers. Together with the broad subharmonic 
peaks observed in the temporal spectrum of Fig. [3J this confirms that the parametric resonance is weakly selective in 
this system. Values of k\ corresponding to significant growth rates are of the same order of magnitude as the primary 
wavenumber kq = 0.82 rad cm" 1 , indicating that the viscosity has a significant effect on the selection of the excited 
resonant triad. For the value of cfq/ f considered in Fig. [6l the maximum growth rate is obtained for the (—,+,—) 
resonance, for &™ ax = 0.75 rad cm -1 . The corresponding predicted wavevector ki is represented as a circle in the 
resonance curve of Fig. [5ja) , and is found in excellent agreement with the experimental measurement of ki (shown 
with an arrow). 

Because of the viscous attenuation, the primary wave amplitude Aq actually depends on the distance from the 
wavemaker. In the measurement area shown in Fig. Eta), spatial variations of ±25% are found around the average 
Ao = 0.29 cm s _1 . Since the growth rate (|30| depends on Ao, this introduces an uncertainty on the predicted value 
of 7, and consequently on the selected secondary wavenumbers. In order to appreciate the influence of the measured 
value of Ao on the predicted triadic resonance, we also plot in Fig. [5] the growth rate of the selected (—,+,—) 
resonance, but for a value of Aq increased by an amount of 25% (continuous line), which corresponds to the wave 
amplitude in the close vicinity of the wavemaker. The maximum growth rate is actually found to strongly depend on 
Aq, with an increase of 30%, indicating that the onset of the parametric instability will take place first close to the 
wavemaker. This strong sensitivity would make any direct comparison with an experimental growth rate too difficult. 
On the other hand, the selected wavenumber £;™ ax is quite robust, showing a slight increase of 6% only when Ao is 
increased by 25%. As a consequence, the uncertainty in the measurement of Ao, which is unavoidable because of the 
viscous attenuation of the primary wave, does not affect significantly the prediction for the most unstable secondary 
wavevectors. 

The size of the circles in Figs. [SJa) and (b) illustrates the uncertainty in the determination of the most unstable 
wavevectors due to the spatial variation of Ao- The relative uncertainty lies in the range 5 — 15% for the range of wave 
frequencies considered here. In spite of this uncertainty, we can conclude that the secondary wavevectors predictions 
from the maximum growth rate criterion are in good agreement with the observed resonant triads. 



C. Dependence of the secondary waves properties on the primary wave frequency 

We finally characterize here the evolution of the secondary wave properties (frequencies and wavenumbers) as the 
frequency of the primary wave is changed. For a given primary wave amplitude Ao, the secondary frequencies ai t 2 
and wavenumbers have been systematically computed according to the maximum growth rate criterion, and 
are reported in Fig. [7] as a function of ctq// G [0,1]. The dotted lines correspond to the (—,—,—) and (—,+,—) 
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FIG. 7. (Color online) Normalized frequencies eri^/oo (a)i an d wavenumbers ^1,2/^0 (b) of the secondary waves, as a function 
of the primary wave frequency (To//. Filled circles and squares with errorbars correspond to experimental measurements. 
Predictions from the triadic resonance instability are represented with dashed thick lines (using absolute maximum growth 
rate criterion) and dotted lines (using maximum growth rate criterion for the ( — ,—,—) and ( — ,+,—) resonances). Predictions 
for the most unstable resonance are ( — , — , — ) for (Jo// < 0.79, and ( — , +, — ) for (Jo// > 0.79. These predictions have been 
computed with a typical amplitude Ao = 0.30 cm s~ : for the primary wave. Continuous solid lines show the allowed range 
around the (— , +, — ) curves, determined by considering an uncertainty of ±50% on Aq. 



resonances, whereas the dashed thick lines are computed from the absolute maximum growth rate among all the 
possible resonances. For (To// > 0.79, the growth rate is maximum on the (— ,±,=p) branch, whereas for ctq// < 0.79 
it is maximum on the (—,—,—) branch. 

In Fig. [71 we also show the experimental measurements of (T1.2 and K1.2 for the range of primary wave frequencies for 
which a subharmonic instability is observed, 0.65 < ao/ f < 0.99. The errorbars show the uncertainties computed from 
the measured frequencies and wavenumbers. The agreement with the predictions from the triadic resonance theory 
is excellent for the (— , +, — ) branch. However, it is not clear why all the measurements actually follow the (— , +, — ) 
branch, although the (— , — , — ) branch is expected to be more unstable for the two data points at ao/ f < 0.79. 

The limited spatial extent of the primary wave along its transverse direction (which represents 4 wavelengths only) 
and its amplitude decay along its propagation direction (because of viscous attenuation) may be responsible for this 
unexpected stability of the (— , — , — ) branch at low (To//. Indeed, the (— , — , — ) branch is associated to wavelengths 
significantly larger than the primary wavelength, so that a large spatial region of nearly homogeneous primary wave 
amplitude is required to sustain such large wavelength secondary waves. On the other hand, the (— , +, — ) resonance 
generates lower wavelengths, which can more easily fit into the limited extent of the primary wave. Finite size effects 
may therefore explain both the preferred (— , +, — ) resonance at er // < 0.79, and the unexpected global stability of 
the primary wave for ao/ f < 0.65. Confinement effects are not described by the present triadic resonance theory, 
which assumes plane waves of infinite spatial extent. Apart from this open issue, we can conclude that, at least for 
sufficiently large forcing frequency, the observed secondary frequencies and wavenumbers are in good quantitative 
agreement with the predictions from the triadic resonance theory. 



V. DISCUSSION AND CONCLUSION 



Using a wavemaker initially designed to generate beams of internal gravity waves in stratified fluids, we have 
successfully generated well-defined plane inertial waves in a rotating water tank. Spectral analysis, performed on 
particle image vclocimetry measurements of this plane inertial wave, has revealed the onset of a parametric instability, 
leading to the emergence of two secondary subharmonic waves. The wavevectors and frequencies of the primary and 
secondary waves are found in good agreement with the spatial and temporal resonance conditions for a resonant triad 
of inertial waves. Moreover, using the triadic resonance theory for inertial waves derived by Smith and Waleffe^ the 
growth rate of the instability has been computed, yielding predictions for the secondary wavevectors and frequencies 
in agreement with the measurements. At low forcing frequency, we observe a departure from these predictions which 
may be associated to the finite size of the primary wave. These finite size effects cannot actually be described within 
the triadic resonant theory, which relies on plane waves of infinite extent. 

Triadic resonant instability for inertial and internal waves share a number of common properties. In particular, 
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equations governing the wave amplitudes equivalent to Eqs. ((27)) and ([28]) may also be derived for a triad of internal 
waves, but in this case, they concern the amplitude of streamfunctions and not of velocities^ The interaction 
coefficients for internal waves G r (with r — 0, 1, 2) can be readily obtained from the interaction coefficients for inertial 
waves C r through a simple exchange of the vertical and horizontal components of the wavevectors, and introducing a 
prefactor: 

C r (k,m) = ^-C r {m,k). (31) 

The KpKqj K r prefactor between the two types of coefficients comes from the fact the wave amplitude is directly given 
by the velocity u in the case of inertial waves, whereas it is given by the strcamfunction ijj ~ u/k in the case of internal 
waves. The exchange of the vertical and horizontal components of the wavevectors comes from the comparison between 
the dispersion relations for inertial and internal waves, cr// = sm/n and a/N = skjn respectively, with f = 10. the 
Coriolis parameter and N the Brunt- Vaisala frequency. The inviscid growth rate of the parametric instability 7 for 
the internal waves is actually equal to the one of inertial waves 7 through 

7 = = ^/ChChKoA = 7, (32) 

where Aq is the primary internal wave amplitude (homogeneous to a streamfunction). Here, the inertial wave amplitude 
Aq (homogeneous to a velocity) identifies with kqAq. This equality between inertial and internal growth rates finally 
shows that the predicted secondary waves should be identical for the two types of waves. 

Interacting inertial waves are of primary importance for the dynamics of rotating turbulence. In the limit of low 
Rossby numbers Ro = U /Q.L, where U and L are characteristic velocity and length scales, rotating turbulence can be 
described as a superposition of weakly interacting inertial waves, whose interactions are directly governed by triadic 
resonances. This is precisely the framework of wave turbulence as analyzed in Refs. HH and|39| in the context of rotating 
turbulence. The parametric instability between three inertial waves can be seen as an elementary process by which 
energy is transferred between wavevectors in rotating turbulence. This anisotropic energy transfer takes place both 
in scales (or wavenumbers) and directions (or angles). The angular energy transfer is always directed towards more 
horizontal wavevectors, providing a clear mechanism by which slow quasi-2D motions become excited^ However, the 
nature of energy transfers through triadic resonance in terms of wavenumbers (or scales) — i.e., whether the energy 
proceeds from large to small scales or inversely — is found to depend on wave amplitude and viscosity. Indeed, it can 
be shown theoretically, within the present triadic resonance framework, that waves of amplitude large compared to ukq 
are unstable with respect to secondary waves of large wavenumbers, producing a direct energy cascade towards small 
scales. On the other hand, waves of amplitude much lower than vkq are found to excite secondary waves of smaller 
wavenumber, hence producing an inverse energy cascade towards larger scales. The net result of this competition is 
delicate to decide, and may contain an answer to the debated issue concerning the direction of the energy cascade in 
rapidly rotating turbulence. 
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